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In this paper we analyze the cycle-to-cycle variations of maximum pressure pmax and 
peak pressure angle apmax in a four-cylinder spark ignition engine. We examine the 
experimental time series of pmax and apmax for three different spark advance angles. 
Using standard statistical techniques such as return maps and histograms we show that 
depending on the spark advance angle, there are significant differences in the fluctuations 
of Pmax and apmax • We also calculate the multiscale entropy of the various time series 
to estimate the effect of randomness in these fluctuations. Finally, we explain how the 
information on both pmax and apmax can be used to develop optimal strategies for 
controlling the combustion process and improving engine performance. 

Keywords: stochastic process; stochastic analysis; nonlinear time series analysis; com- 
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1. Introduction 

Instabilities in combustion processes were observed from the very beginning of the 
Spark Ignition (SI) engine development These instabilities may cause fluctua- 
tions in the power output as well as fluctuations of the burned fuel mass. As a 
consequence, the mean effective torque may be reduced by as much as 20%. It 
is therefore not surprising that the instabilities were identified as a fundamental 
combustion problem in spark ignition engines . A disturbing feature of these in- 
stabilities is the u npredictab le character of their occurrence, which makes an engine 
difficult to control I Sl^ l ^ l ^ l ^ i Recognition and elimination of their so urces hav e been 
one of the main issues in SI engine technology in the last century [ 8 | 9 | 1Q | 11 ] How- 
ever, in spite of great efforts made in clarifying the various aspects of combustion 
instabilities in the past several years, the proble m of a stable combustion process 
control in an SI engine has not yet been The main factors of com- 

bustion instabilities as classified by Heywood ^ are aerodynamics in the cylinder 
during combustion, amounts of fuel, air and recycled exhaust gases supplied to the 
cylinder, and composition of the local mixture near the spark plug. Among the var- 
ious researchers, Winsor et a/. explored the turbulent aspects of combustion and 
studied the nature of pressure fluctuations in an SI engine. Recent attempts have 
focused on the development of nonlinear dynamical models of the combustion pro- 
cess. For example. Daily ^1^1 fQ^j^iulated a simple nonlinear model and demonstrated 
that pressure variations of a chaotic type can originate from the initial conditions 
prevailing at the beginning of each cycle. Using a nonlinear model, Kantor ^ an- 
alyzed the cycle-to-cycle variations of the process variables including combustion 
temperature. Foakes and Pollard , Chew et a/, ^^and Letellier et a/.H^ examined 
nonlinear models of pressure variations and calculated their Lyapunov exponents 

Subsequently, Daw et al included the effect of exhaust gas circulation in a 

nonlinear model to investigate engine dynamics. But the most convincing argument 
for the presence of nonlinear dynamics in the combustion process came from the 
concept of time irreversibility of heat release as shown by Green et al. , Wagner 
et al. ttOl and Daw et al. More recently, on the basis of their analysis of cycle-to- 
cycle pressure variations, Wendeker et al. proposed an intermittency mechanism 
leading to chaotic combustion. 

In this paper we continue these studies with an examination of the experimental 
time series of maximum pressure Pmax and peak pressure angle apmax for three 
different spark advance angles. Our paper is divided into four sections. After this 
introductory section (Sec. 1), we provide a description of our experimental standing 
in Sec. 2 and outline the procedure of performing measurements. This is followed 
by an analysis of the statistical properties of the time series of Pmax and apmax 
using return maps and histograms. We also calculate the multiscale entropy (MSE) 
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Fig. 1. Schematic diagram of the experimental stand for pressure measurements. 



of the various time series in order to estimate their complexity and dehneate the 
effect of randomness in the fluctuations of Pmax and apmax- Finally, in Sec. 4, we 
discuss the results of our analysis and use them to propose a criterion for optimal 
efficiency of combustion. 



2. Experimental Stand for Pressure Measurement 

In our experimental stand (Fig. [T]), the pressure inside the cylinder is measured by 
a standard method using a piezo-electric sensor. It should be noted that pressure is 
the best known quantity that can be directly measured to analyze engine dynamics 
| 5 | 8 > 26 27 28 i 29l Actual cylinder pressure together with the internal cylinder volume 
can be used to obtain the indicated mean effective pressure (IMEP), to calculate 
the engine torque and efficiency, and also to estimate the magnitudes of important 
process variables such as burning rate, bulk temperature and heat release. Further- 
more, statistical analysis of internal pressure data can provide information about 
the stability of the combustion process. 

Our equipment provides a direct measure of the pressure inside the cylinder of 
a spark ignition engine. Internal pressure data were recorded at the Engine Labora- 
tory of the Technical University of Lublin, where we conducted a series of tests. The 
pressure traces were generated on a 1998 cm^ Holden 2.0 MPFI engine running at 
1000 revolutions per minute. Measurements were made with a resolution of about 
0.7^ of crankshaft revolution. A single experiment consisted of about 2000 cycles 
of engine work. Data were collected at different spark timings with spark advance 
angles of 5, 15 and 30 degrees before the top dead center (TDC). The engine speed, 
fuel-air ratio, and throttle setting were all held constant throughout the data col- 
lection period. Intake air pressure was maintained at 40 kPa and a stoichiometric 
mixture was used. 

A few examples of the pressure measurements are presented in Fig. O Here we 
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Fig. 2. Variations of internal pressure with crank angle during the first 20 combustion cycles. The 
panels (a), (b) and (c) correspond to the spark advance angles Aaz = 5°, 15° and 30°, respectively. 



have plotted the variation of pressure with crank angle during the first 20 com- 
bustion cycles for three different values of the spark advance angles: Aaz = 5^, 
15^ and 30^. Note that the pressure measurements reflect the combined effects of 
cylinder volume compression and combustion of the fuel-air mixture. From these 
measurements the values of maximum pressure Pmax and the peak pressure angle 
(^pmax can be easily identified. 

Before we go into a detailed analysis of the time series of Pmax and apmax in 
consecutive combustion cycles, we would like to point out the importance of these 
quantities for combustion control and diagnostics. Note that the output torque 

usually scales as the square of the mean values of Pmax or apmax for different 
speeds 
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Fig. 3. Time series of cycle-to-cycle variations of maximum pressure pmax (in black) and apmax 
(in grey) for three values of the spark advance angle: Aaz = 5°, 15° and 30°, for Figs. [3ji-c, 
respectively. 



3. Cyclic Time Series and Pressure Fluctuations 

From Fig. [2] we observe that the cycle-to-cycle fluctuations of pressure are present in 
all cases considered and they change with the spark advance angle Aa^. It should 
be noted that the magnitude of peak pressure, which is directly related to the 
power output, depends on Aa^; Pmax is higher for larger Aa^. Unfortunately, the 
advantage of maximum pressure increase is offset by the disadvantage of increasing 
fluctuations. In terms of the peak pressure Pmax^ this effect can be observed in Figs. 
[3^-c. Here we have also plotted the time series of the corresponding angles apmax 
for which the peak pressure Pmax is reached. Surprisingly, changes in combustion 
conditions which lead to increased fluctuations of Pmax are associated with decreased 
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fluctuations of a^. 



^pmax • 



3.1. Return Maps 

To explore the properties of the time series of Pmax and apmax in niore detail, we 
have plotted their return maps for the various spark advance angles. These are 
shown in Fig. [H We define a sequential transformation for any cycle i: 

+ (3.1) 

where F can be identified with Pmax or apmax- In a return map, black points in 

the coordinates {F{i-\- l),F(z)) are marked. Such maps are very useful in examining 

changes of the corresponding quantity F{i). The method of return maps has become 

a standard tool in the investigation of engine dynamics but they are usually applied 

90 99 ^0 

to burned mass ratio or heat release per cycle . Let us analyze the first dia- 

gram (Fig. 13^) which corresponds to the pressure history in Fig. [2^ (Aa^ = 5^). The 
characteristic black right angle in the left bottom corner of the return map indicates 
a very low level of combustion or misfire instabilities. The corner point represents 
the lowest value of Pmax which can be identified as peak compression pressure. In 
terms of peak pressure angles apmax (Fig- Sb) we observe that the points lying at 
the lowest angular position around 360^ (for Aa^ = 5^) correspond to a maximum 
compression point - the Top Dead Point (TPD). However, in this case (Aa^ = S*^) 
we cannot draw any conclusion about misfires because w e are not an alyzing heat 
release here. Extended studies of this effect can be found inl22ll0l30l32]^ Figs. [It, d 
and Figs. [4t,f, we show the return maps for the spark advance angles Aa^ = 15^ 
and 30^, respectively. In these cases the fluctuations may have a different origin. 
Following the previous analysis for Aa^ = 15^ (Figs. [It, d), we have detected a sin- 
gular point with relatively weak combustion, while for Aa^ = 30^ such points have 
disappeared. This means that the identification of pressure peak caused by combus- 
tion is certain in that regime. Of course, the measured fluctuations will change with 
changes in the spark advance angle Aaz Statistical features covered by the re- 
turn maps are in agreement with the time series histories depicted in Fig.[3ti-c. Note 
that the fluctuations in Pmax increase while those in apmax decrease with increasing 
Aaz. We should also mention that all the maps possesses the approximate diago- 
nal {F{i + 1) versus F{i)) symmetry; this can be interpreted as the time reversal 
symmetry of the corresponding time series of Pmax and apmax- This symmetry has 
not been broken in contrast to other reports on combustion experiments (Green at 
al. 1221) . The main difference is that Green at al. 1221 examined the lean combustion 
limit whereas our fuel- air mixture is stoichiometric. Clearly, preserved time-reversal 
symmetry may appear if the stochastic component is significant. Alternatively, for 
deterministic but chaotic time series one usually observes asymmetry in the return 
maps '22i associated with a broken time-reversal symmetry. 
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Fig. 4. Return maps of cyclic pressure maxima pmax and the corresponding angles apmax for three 
spark advance angles Aaz denoted by SA= 5°, 15° and 30°. 
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Fig. 5. Histograms of cyclic pressure maxima Pmax and the corresponding angles OLpmax 
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3.2. Histograms 

Here we perform further analysis of the Prnax and apmax time series using his- 
tograms. From the time series we calculate the probability P^- 

= 7V ^ " - ^n+l), (3.2) 
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where Q{x) is the Heaviside step function, F{i) is a sequentially measured quantity, 
and Fn denotes a discretized value. Here N total = 1991 is the total number of 
combustion cycles in a single time series. In Figs. [5^-1 we have plotted the histograms 
of the time series for Pmax or apmax for the various cases considered in Figs. [4^- 
f. Interestingly, the combustion process represented by the histograms in Figs. [5^ 
and b can be identified as a process with an ambiguous position of maximum. 
This effect was also found in the paper by Litak et al. In fact the double peak 
structure is clearly visible in Fig. [Sb- They are directly related to compression and 
combustion phenomena. Note that the exponential distribution in Fig. [5^1 indicates 
that the measurement process is not continuous, similar to the situation in a cutting 
process with additional cracking effects In our case, identification of the pressure 
peak is not certain because from time to time, especially for weak combustion it is 
interchanged with compression. On the basis of our experimental data we cannot 
entirely exclude the possibility of misfires In Figs.[5fc-d and Figs.[5fe-f, we show 
the distributions of Pmax and apmax for Aa^ = 15^ and Aa^ = 30^, respectively. 
Comparing these figures with Figs.[5K-b, we see that they are substantially different 
having more or less the Gaussian shape with small asymmetry. Again we observe a 
characteristic broadening in Pmax and a narrowing in aprno^with increasing Aa^ . 

3.3. Multiscale Entropy 

Many physical systems evolve on multiple temporal and/or spatial scales and are 
governed by complex dynamics. In recent years there has been a great deal of in- 
terest in quantifying the complexity of these systems. However, no clear and unam- 
biguous definition of complexity has been established in the literature. Intuitively, 

OA 

complexity is associated with "meaningful structural richness" Several notions 
of entropy have been introduced to describe complexity in a more precise way. But, 
as Costa et al.^^have pointed out, these traditional entropy measures quantify the 
regularity (predictability) of the underlying time series on a single scale, and there 
is no direct correspondence between regularity and complexity. For instance, neither 
completely predictable (e.g., periodic) signals, which have minimum entropy, nor 
completely unpredictable (e.g., uncorrelated random) signals, which have maximum 
entropy, are truly complex. For complex systems with multiple temporal or spatial 
scales, a definition of complexity should include these multiscale features. Recently 
Costa et al. introduced the concept of multiscale entropy (MSE) to describe the 
complexity of such multiscale systems. They have successfully used this concept to 
describe the nature of complexity in physiological time series such as those associ- 
ated with cardiac dynamics ^3 and gait mechanics ESl. In this paper we perform a 
multiscale entropy analysis as a measure of complexity in cycle-to-cycle variations 
of maximum pressure and peak pressure angle in a spark ignition engine. MSE is 
based on a coarse-graining procedure and can be carried out on a time series as 
follows. 

For a given time series {Fi, F2, F^}, where = F(i), multiple coarse-grained 
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Fig. 6. Schematic diagram of the coarse-graining procedure used in the calculation of multiscale 
sample entropy (after [34]). Note, F^'^'^ = Fi and r is the scale factor. 




Fig. 7. Sample entropy Se calculated for normalized time series Fi = pmax{i) / cf{Pmax) (a) and 
apmax / o-{apmax) (b) versus the scale factor r. (j(pmax) and cr(apmax) denote square root devi- 
ations of Pmax and apmax, respectively. The length of chain m = 3 and similarity factor r = 0.15 
for all time series. Curves '1', '2' and '3' correspond to Aaz = 5^, 15^ and 30^, respectively. 



time series are constructed by averaging the data points within non-overlapping 
windows of increasing length as shown in Fig. [6l Each element of the coarse-grained 
time series is computed according to the equation: 



1 



(3.3) 



i=(j-l)r + l 



Here r represents the scale factor, and 1 < j < N/r^ the length of each coarse- 
grained time series being equal to N/r. Note that for r = 1, the coarse-grained time 
series is simply the original time series. 

The sample entropy Se^ proposed originally by Richman and Moorman ^^H] 
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calculated for the sequence of m consecutive data points which are 'similar' to each 
other and will remain similar when one or more consecutive data points are included 
ESI. Here 'similar' means that the value of a specific measure of distance is less than 
a prescribed amount r. Therefore, the sample entropy depends on two parameters 
m and r. We have 

SE{m, r, N) = In (3.4) 

The sample entropy has been calculated for each of the coarse-grained time series 
of Pmax and apmax ^ud then plotted as a function of the scale factor (r). The results 
are shown in Fig. [71 where Fig. [7^ corresponds to Pmax and Fig. [Tb to apmax- It 
should be noted that for a simple but stochastic system the sample entropy Se has a 
hyperbolic shape. This shape can be identified only for curve '2' in Fig.[7K while the 
curves '1' and '3' are definitely more fiat. In Fig. [TJd, on the other hand, all curves 
have a mild dependence on r, and can be described approximately by horizontal 
lines. Such horizontal lines are typical for strongly asymmetric distributions such 
as 1/f noise. In fact, if one looks more carefully at the slopes of the histograms of 
our signals, one can notice the obvious asymmetry in Figs. and e, and to some 
extent in Figs.[5]i,f. Note that the peculiar nonlinear dynamics and a limited length 
of time series can additionally disturb these characteristics and produce fiuctuations 
along the r axis (Fig. [7j). Interestingly the results show that Pmax and apmax are 
characterized by quite different aspects of the same dynamical process. 



4. Discussion and Conclusions 

We have used the maximum pressure Pmax and the peak pressure angle apmax 
as variables for describing the combustion process in a spark ignition engine. Their 
cycle-to-cycle variations are analyzed using different techniques such as return maps, 
histograms and multiscale entropy. Until now, researchers have used only one of 
these variables to formulate an objective function for optimization and control of 
engine dynamics. Prompte d by the se complementary analyses of Pmax and apmax 
presented in the literature l ^^ l ^^ l ^^ l^ we propose to use a mutual criterion for com- 
bustion efficiency in the SI engine. Instead of minimizing the fluctuations of Pmax 
or apmax by looking at their smallest standard deviation cr{pmax)^ or cr (apmax)- 



(a) cr{pmax) = 

(b) cr{apmax) 



^Pmax \ 
1 

Sanmnx 



Ntotal 

IT. V {Pmax{i) 

^^total 



\2 

Pmax ) 



(4.1) 



-j^ ' total 

\ ^ ^ i^pmaxi,'^) apmax) 
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Fig. 8. Square deviations crpmax (a), (j(apmax) (b) and cr(pmax, ctpmax) (c) as a generalized 
efficiency criterion (notation (a)-(c) as in Eqs. I4.H4.2|) . 



we look for a two- variable square deviation (j{pmax^ Oipmax) given by. 



\ 




(i)-Pmax) , {(^pmax{i) 



^total 



(4.2) 



from all 



where Spmax and Sapmax are chosen as the largest of (ipmax and (T{apr 
examined cases and are treated as normalization. 

All the above square deviations are plotted in Fig. [H In this figure note that 
the plot for (j{pmax^ crpmax) has a minimum when Aaz = 15^. In our previous 
papeiflSl^ where we analyzed noise effects in the cycle-to-cycle fiuctuations of heat 
release, we noticed that Aaz = 15^ is the case of optimal combustion conditions. 
Now we see that for Aaz = 30^ the oscillations of Pmax are the largest. On the 
other hand, for Aaz = 5^ the fluctuations in cr{apmax) are significant. The average 
heat release which varies with the spark advance angle has magnitudes of 352. 54 J 
for Aaz = 5^ 396.67J for Aaz = 15^ and 377.04J for Aaz = 30^ It is the 
largest for Aaz = 15^ indicating the largest burning rate of fuel while the output 
torque, for the same speed of the crankshaft, is changing from the largest value: 
5* = 30Nm, in the case of Aaz = 30^, to a slightly smaller value: 5* = 28Nm for 
Aaz = 15^, and to an even smaller value: S = 21Nm for Aaz = 5^. Knowing 
that the fresh fuel rate was the same in all the cases, we concluded that there were 
better combustion conditions for larger spark advance angles. Additionally, in case 
of Aaz = 5^ some amount of fuel was not burned because of weak combustion or 
misfires. In the previous paper our analysis of noise level led us to the conclusion 
that in the case of Aaz = 15^, the combustion process was most stable. 

Here we have found another argument leading to the same conclusion based on 
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Fig. 9. Maximum pressure-angle diagrams (pmax, c^pmax) for three advance angles (SA = Aaz 
5° (a), 15° (b), 30° (c)). 



the two variables Pmax and apmax- Furthermore, from the multiscale entropy analy- 
sis we find that the Pmax sequence is governed by the simplest dynamics. Examining 
both [pmax^ Oipmax] time scrics simultaneously it becomes clear that the two variable 
square deviation {cF{praax-,ocpmax)) may play an important role in the optimization 
procedure. In Figs.[9ti-c we plotted the variables for each combustion cycle. Follow- 
ing our philosophy, these plots may be called 'a bifurcation diagram'. They show 
the transition from the case of Ofp^aa^-dominated fluctuations (Fig. [9^1 for Aa^ = 5^) 
to the case of prna^c-dominated fluctuations (Fig. [9b for Aa^ = 30^) through the 
intermediate case (Fig.[9]3 for Aa^ = 15^). One can easily see the qualitative change 
in the fluctuation type starting from a vertical line structure in Fig. [9^ and ending 
with an almost horizontal line in Fig. [9b. It is worth pointing out that both types of 
fluctuations are harmful for a stable engine operation and may result in decreased 
power output. It appears that the intermediate case is a compromise between these 
two types of fluctuations. The detailed analysis of the relations between Prnax and 
o^pmax can provide further useful information on combustion dynamics as the rate 
of slow and/or fast individual burning processes during engine use corresponding 
to slow and fast heat release. The two extreme realizations of a combustion pro- 
cess presented in Fig. [9^1 and Fig. [9fc (for Aa^ = 5^ and 30^, respectively) can be 
directly related to slow and fast burning Obviously the fast burning is not so 
effective as slow burning because it may involve some events corresponding to a 
spatially limited engine volume. Thus, as expected, if slow burning dominates, then 
the fluctuations of Pmax are small (Figs. [9^1 and[3ti), whereas for fast burning the 
fluctuations of Pmax are larger (Figs. [9t and[3t). Interestingly, measured fluctua- 
tions of aprnax are governed mainly by the ambiguity of identiflcation of maximum 
pressure and possibly by a misflre phenomenon These two effects have oppo- 
site influences on Pmax and a^rnax- On this basis our formulation of optimization 
procedure (Eq. 14. 2p seems to be a natural consequence of a detailed analysis of the 
pressure-angle diagrams Figs. [9K-c. The minimum standard deviation deflned in 
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the {pmax: c^pmax) plane, realized for Aaz = 15^, can be regarded as a compromise 
between the maximizing of power output and minimizing of pressure fluctuations. 
Clearly, larger power is accompanied by unstable combustion with larger pressure 
oscillations. 

It is worth noting that our simple mutual criterion seems to unify these ap- 
proaches. Using this criterion we can associate the optimal combustion conditions 
with the spark advance angle Aaz = 15^. Of course, for any practical application 
we would need to consider many more values of Aa^ and more importantly check 
the output torque^. 
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